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In this paper, we study theoretically the bipartite and tripartite continuous variable entanglement 
as well as the normal-mode splitting in a single-atom cavity optomechanical system with intensity- 
dependent coupling. The system under consideration is formed by a Fabry-Perot cavity with a thin 
vibrating end mirror and a two-level atom in the Gaussian standing-wave of the cavity mode. We 
first derive the general form of Hamiltonian describing the tripartite intensity-dependent atom-field- 
mirror coupling due to the presence of cavity mode structure. We then restrict our treatment to the 
first vibrational sideband of the mechanical resonator and derive a novel form of tripartite atom- 
field-mirror Hamiltonian. We show that when the optical cavity is intensely driven one can generate 
bipartite entanglement between any pair of the tripartite system, and that, due to entanglement 
sharing, the atom-mirror entanglement is efficiently generated at the expense of optical-mechanical 
and optical-atom entanglement. We also find that in such a system, when the Lamb-Dicke parameter 
is large enough one can simultaneously observe the normal mode splitting into three modes. 

PACS numbers: 37.30.+i, 03.67.Bg, 42.50.Wk, 85.85.+J 



I. INTRODUCTION 

Cavity optomechanics is a rapidly growing field of re- 
search that is concerned with the interaction between 
a mechanical resonator (MR) and the radiation pres- 
sure of an optical cavity field [THE]. The optomechan- 
ical coupling widely employed for a large variety of 
applications 7 , more commonly as a sensor for the de- 
tection of weak forces [H] and small displacements [3] or 
an actuator in integrated electrical, optical, and opto- 
electronical systems [101 111) . However, the most exper- 
imental and theoretical efforts are devoted to cooling 
and trapping such mechanical resonators to their quan- 
tum ground state, which more recently have been done 
successfully [T2. Furthermore, in Ref. [13] the authors 
have proposed a different scheme to enhance the cooling 
process by using the photothermal (bolometric) force [H]. 
They have taken into account the noise effects due to the 
granular nature of photon absorption and finally have 
shown that the mechanical resonator can achieve the low- 
est phonon occupation number by means of this proce- 
dure. Moreover, it seems promising for the realization of 
long-range interaction between qubits in future quantum 
information hardwares |15) , and for probing quantum me- 
chanics at increasingly large mass and length scales [IB]. 
The coupling of a MR via radiation pressure to a cavity 
field shows interesting similarities to an intracavity non- 
linear Kerr-like interaction [5] or even a more complicated 
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form of nonlinearitv|17j. 

To observe and control quantum behavior in an op- 
tomechanical system, it is essential to increase the 
strength of the coupling between the mechanical and 
optical degree of freedom. However, the form of this 
coupling(e.g., linear or nonlinear) is crucial in determin- 
ing which phenomena can be observed in such a sys- 
tem. Thanks to rapid progress of nano-tcchnology, it 
has been possible to manipulate the optomechanical cou- 
pling in quantum optomechanical hybrid systems. In this 
direction most experimental and theoretical efforts are 
devoted to entangling a MR either with a single atom 
[TSH2"2"] or with atomic ensembles [2"3Tf2"5] . entangling a 
nanomechanical oscillator with a Cooper-pair box [29] . 
and entangling two charge qubits [30] or two Josephson 
junctions 131] via nanomechanical resonators. Alterna- 
tively, schemes for entangling a superconducting coplanar 
waveguide field with a nanomechanical resonator, either 
via a Cooper-pair box within the waveguide [32] , or via 
direct capacitive coupling [33], have been proposed. 

In Ref. [23] the authors have proposed a scheme for the 
realization of a hybrid, strongly quantum-correlated sys- 
tem consisting of an atomic ensemble surrounded by a 
high-finesse optical cavity with a vibrating mirror. They 
have shown that, in an experimentally accessible pa- 
rameter regime, the steady state of the system shows 
both tripartite and bipartite continuous variable(CV) en- 
tanglement. More recently, the dynamics of a movable 
mirror of a cavity coupled through radiation pressure 
to the light scattered from ultracold atoms in an opti- 
cal lattice has been investigated Ref .[34 . The author 
has shown that in the presence of atom-atom interaction 
as a source of nonlinearitv|35j. the coupling of the me- 
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chanical oscillator, the cavity field fluctuations and the 
condensate fluctuations (Bogoliubov mode) leads to the 
splitting of the normal mode into three modes (normal- 
mode splitting(NMS) 36-40]). The system described 
there shows a complex interplay between distinctly three 
systems namely, the nanomechanical cantilever, optical 
microcavity and the gas of ultracold atoms. 

The optomechanical NMS is one of the fascinating phe- 
nomena arising from the strong coupling between the cav- 
ity and the mechanical mirror [4 1 -43J. In Rcf.[42] it has 
been shown that the cooling of mechanical oscillators in 
the resolved sideband regime at high driving power laser 
can entail the appearance of NMS. Moreover, the dy- 
namics of a movable mirror of a nonlinear optical cavity 
is considered in Ref . |43j . It has been shown that a \^ 
medium with a strong Kerr nonlinearity placed inside 
the cavity inhibits the NMS due to the photon blockade 
mechanism(this just happens only if the Kerr nonlinear- 
ity is much greater than the cavity decay rate). As the 
authors have shown in Refs.[Mj and [13] the nonlinear- 
ity plays a crucial role in the appearance of NMS in the 
optomechnical systems. 

The main purpose of the present paper is to study 
the quantum behavior of an atom-assisted cavity op- 
tomechanical system in which a single two-level atom is 
trapped in the standing-wave light field of a single-port 
Fabry-Perot cavity. The infinite set of optical modes of 
the cavity can be described by Hermit-Gauss modes. As 
we will see, the intracavity mode structure can be em- 
ployed to realize a type of intensity-dependent coupling 
of the single atom to the vibrational mode of MR. The 
presence of such intensity-dependent interaction modifies 
the dynamics of the system, the entanglement properties 
and the displacement spectrum of MR. We show that in 
the first vibrational sideband of MR, the stationary, i.e., 
long-lived, atom-mirror entanglement can be generated 
by proper matching the Lamb-Dicke parameter(LDP). 
This parameter plays an important role in our investi- 
gation in the sense that it determines the strength of 
the nonlinearity in the system. We show that the bi- 
partite entanglement between the subsystems extremely 
depends on the LDP. It is also remarkable that, in the 
steady-state condition, the high resolution of NMS in the 
form of three-mode splitting is approached. In particu- 
lar, the appearance of intensity-dependent coupling leads 
to a progressive increasing of NMS due to the strong non- 
linear atom-field-mirror interaction. 

The paper is organized as follows. In Sec. II we derive 
an intensity-dependent Hamiltonian describing the triple 
coupling of atom-field-mirror through j-phonon excita- 
tions of the vibrational sideband. In Sec. Ill, we derive 
the quantum Langevin equations (QLEs) and linearize 
them around the semiclassical steady state. In Sec. IV 
we study the steady state of the system and quantify 
the entanglement properties of the system by using the 
logarithmic negativity. In Sec. V we investigate the ap- 
pearance of NMS in the displacement spectrum of the 
mirror. Our conclusions are summarised in Sec. VI. 



II. MODEL 

The system studied in this paper is sketched in Fig. 1. 
It consists of a hybrid system formed by a single two- 
level atom with transition frequency lo c which trapped in 
the standing-wave light field of a single-port Fabry-Perot 
cavity with a movable mirror coated on the plane side of 
a mechanical resonator. The geometry of the resonator 
determines the spatial structure of the acoustic modes. 
The movable mirror is treated as a quantum mechani- 
cal harmonic oscillator with effective mass m, frequency 
oj m , and energy decay rate j m . The system is also coher- 
ently driven by a laser field with frequency uii through 
the cavity mirror with amplitude £ . We assume that the 
single atom is indirectly coupled to the mechanical os- 
cillator via the common interaction with the intracavity 
field with frequency lu c . 



Optical cavity 




da 



FIG. 1. The schematic of the atom-assisted optomechanical 
system. It contains an optical cavity ended with a fixed mirror 
and a slightly moving mirror which is attached to a spring. 
Inside the cavity there is a two-level atom. The system is 
coherently driven by a laser field 

In our investigation we can restrict the model to the 
case of single-cavity and mechanical modes. This is jus- 
tified when the cavity free spectral range is much larger 
than the mechanical frequency <j m (i.e., not too large cav- 
ities). In this case, scattering of photons from the driven 
mode into other cavity modes is negligible [44 and the in- 
put laser successfully drives only one cavity mode. This 
guarantees the fact that only one cavity mode partici- 
pates in the optomechanical interaction and the neigh- 
bouring modes are not excited by a single central fre- 
quency input laser. In addition, one can restrict to a 
single mechanical mode when the detection bandwidth is 
chosen such that it includes only a single, isolated, me- 
chanical resonance and the mode-mode coupling is neg- 
ligible gS]. 
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A. Hamiltonian of the system 

In the absence of dissipation and fluctuations, the total 
Hamiltonian of the system is given by the sum of three 
term; the free evolution termfini HI] 

Hq = hu] c a) a + hm^b^b H 2~ crZ ' W 

the interaction term 



-h£ (b + tf)^ a + hx m ni(ro, x)[aa + +h.c], (2) 



and the laser driven term 



iJdri = ih£(a'e 



), 



(3) 



where a([a, at] = 1) is the annihilation operator of the 
cavity field with the decay rate k, b([b, &t] = 1) is the 
motional annihilation operator of the MR, and the single 
two-level atom is described by the spin-1/2 algebra of the 
Pauli matrices a~ , <j + and <r z which satisfy the commu- 
tation relations [a + , a~] = <r z and [er^er*] = ±2cr . It 
should be noted that the free HamiltonianQ has been 
written within the Raman-Nath approximation [5T] . i.e., 
in the limit when the atom is allowed only to move over 
a distance which is much less than wavelength of the 
light. Therefore, in this approximation, one can ne- 
glect the kinetic energy of the atom. The first term of 
flint is the optomechanical coupling with the radiation- 
pressure coupling constant £ = ( w c/ 'L)xzp-p, m which 
xzpf = \J h/moj m is the zero point fluctuations of me- 
chanical oscillator. The second term of H lnt denotes the 
"three-body" interactions among the atom, the cavity 
field and the vibration of the mirror. The field-atom 
coupling rate in terms of an infinite set of optical modes 
is well described by the Hermite-Gauss modes [151 |4"T] 
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where, for m, n = 0, 1..., I = 1, 2, 
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w{x)Vir2 n+m - 2 mln\L 



, (4) 



, (5) 



z + y 

ipmni(x,y,z) = kx - <j)(x)(m + n + 1) + fc ,(6) 

2R[x) 



Here go = uJ^/eoV, €q is the vacuum permittivity, V 
shows the volume of the cavity and fj, is the electric-dipole 
transition matrix clement. H n (y) is the n-th Hermite 



polynomial, w(x) = w 1 + (f-) 



is the beam waist 



at x which is defined as the distance out from the axis 
center of the beam where the irradiance drops to 1/e 2 
of its values on axis, R(x) — x + x 2 R /x is the radius of 
curvature of the wavefront at x, 4>{x) — aicta.n(x/xR) 
is the Gouy phase shift 146], wq is the cavity waist ra- 
dius which depends on the geometry of the Fabry-Perot 
cavity and xr = uigfc/2 is the Rayleigh range which 



combines the wavelength and waist radius into a sin- 
gle parameter and completely describes the divergence 
of the Gaussian beam. Note that the Rayleigh range 
is the distance from the beam waist to the point at 
which the beam radius has increased to y/2wo- The 
coupling rate Xmni(ro,x) depends on the initial atomic 
position r (measured from the cavity waist) as well as 
the displacement x = xzpf(& + & ) of the mirror due to 
k = cu c ff(x)/c, where ui g(x) = w c (l — j). As we will 
see in the next section, this dependence to the position 
of MR is responsible for the appearance of a new type 
of optomechanical nonlinearity. Finally, the Hamiltonian 
-ffdri describes the input driving by a laser with frequency 
uji and amplitude \£ \ = 2kP / huii , where P is the input 
laser power and k is the cavity loss rate through its input 
port. 



B. The nonlinear atom-field-mirror coupling 

As we have seen the Gaussian standing-wave struc- 
ture of the cavity mode leads to the field-atom coupling 
rate Xmniifo, x). Such field-atom coupling in the pres- 
ence of the mode structure of the field has been stud- 
ied extensively in the literature and it has been shown 
that a certain type of nonlinearity is prepared in the 
field-atom system. For instance, in Refs. [4"%1 [4"5] the in- 
fluences of the atomic motion and the field-mode struc- 
ture on the atomic dynamics have been investigated. It 
has been shown that the atomic motion and the field- 
mode structure give rise to nonlinear transient effects in 
the atomic population which are similar to self-induced 
transparency and adiabatic effects. In our treatment, 
the spatial field-mode structure leads to the appearance 
of an intensity-dependent interaction among the intra- 
cavity optical mode, the MR and the single atom. To 
show this, we assume that the atom is well located at the 
transverse (polar) coordinate(measured from the cylin- 
drically symmetric cavity axis along the x direction) 



Pa = V^o + yo = /J>w(xo) where < fi < 1. In the 
x direction the localization of the atom can be expressed 
as fco^o — £7r for e > 0, where fco = w c /c- 

At lowest order of the optical modes, i.e., m = n — 
0, 1 = 1, the tripartite coupling rate reduces to 



Xqo\(?q,x) = x{x) 



2.9o 



e fJ, w(xo)v / TrL 



7T 2(IX - 

" Mxo) -2 + J^. 



. . (7) 

which can be rewritten in terms of the mirror position 
by using the position dependence of the wavelength k = 
fc (l — x/L) as 



Xi x ) = / — sm 

e fJ "w(xo)V'n'L 



TJ X 



where % 



W„ fen L 



and 



= (1 



2//, 



2 , 2)^0 - 4>(X ) - 
0^0 



(8) 



(9) 
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By substituting x — xzpF(b + b^) in Eq.Q we obtain 
(&,6 f ) = / g ° N ^_ (e' e exp[»? 7 (6 + b t )]-/i.c} ; (10) 



ie tJ, w(xo)y / VrZ 
where the parameter 



X 



w 2 k 2 L 



mui m 



(11) 



is the so-called "Lamb-Dicke parameter" (LDP). By using 
the Baker-Campbell-Hausdorff theorem in Eq.(10) and 
expanding the exponential terms in terms of b and w, 
the coupling rate can be written as 



X 



5oe 



V/2 



ttL y 



jO 



{irjtf)™' (irjby 
m\m!\ 



-h.cX. 
(12) 

By using the bosonic commutation relation of the oper- 
ators b and w , the j-th term of the field-atom coupling 
rate is obtained as follows 

Xj(b,V)= 



' m).{m + jV. 



- /l. 



where o 7 - u = g ° e . — ^!?L e' e describes the effective atom- 

field-mirror coupling rate and the Hermitian nonlinearity 
function fj{rib) defined by 



/iK) 



(irf) 2m n b \ 



»b- 



m\(m + j)\(nb — m)\ (*n.6 + i)! 



(14) 

with rib — b^b and L J nb (—i] 2 ) as the associated Laguerre 
polynomial, describes a nonlinear atom-field-mirror cou- 
pling through j-phonon excitations of the vibrational 
sideband. The nonlinearity function fj{rib) has a cen- 
tral role in our treatment. It determines the form of 
nonlinearity of the intensity dependence of the coupling 
among the cavity field, the MR and the single atom. As 
we will see, this function drastically influences the dy- 
namics of the system, its entanglement properties and 
its responsible for the appearance of NMS with high vis- 
ibility in the displacement spectrum of the MR. Fig. 2(a) 
shows the nonlinearity function fj(rib) as a function of 
n b and for different values of phonon excitation number 
j. As is seen, this function has maximum contribution 
around small values of vibrational acoustic excitation rib- 
Furthermore, by increasing the number j the strength of 
the nonlinearity function fj(rib) decreases considerably. 
On the other hand, Fig. 2(b) shows that the nonlinearity 
decreases by increasing the LDP. It is remarkable that, 
for the higher orders of the vibrational sideband, j > 3, 
this function is not any more sensitive to the LDP. 

Now, by substituting Eq.(13) into the interaction 
Hamiltonian of Eq. ([2]), we obtain the nonlinear form 
of the Hamiltonian as 



H 



U) 



h^ (b+b r )a 1i a+h s J >(6 t ) j fj(n b )+h.c 



+h.c 
(15)" 
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FIG. 2. The nonlinearity function /j(ri(,) as a function of:(a) 
phonon number rib for r\ = 0.08 and for different values of 
vibrational sideband j and (b) the Lamb-Dicke parameter, 77, 
for rib = 10 and for different values for vibrational sideband 
j- 



Near the photon-phonon resonance |19| where the fre- 
quencies satisfy u m + uj c — co c ~ 0, the rotating- wave 
approximation reduces the above Hamiltonian to 



H£ (b + tf)a? a + H gj^(tf) j j j (n b )a(T + + h.c 



(16) 

This Hamiltonian describes a nonlinear tripartite atom- 
field-mirror coupling and represents a novel type of 
optomechanical intensity-dependent interaction. The 
Hamiltonian ( j 16[ ) is general and one can recover the re- 
sults of Ref.|24j by taking j = 0, fj{rib) — > 1 or results 
of Refs.pl] and [21] by setting j = 1, f 3 {nb) 1. To 
study the system dynamics we restrict our investigation 
by considering the first excitation of the vibrational side- 
band i.e., j = 1. In this limit one may use the following 
simple form of the Hamiltonian 



0=i) 



-K£o(b+b ')a 'a+hgfj, b'f(rib)aa 



-a a} f(n b )b 
(17) 
Since we 



where f{n b ) = fi(n b ) and g^ = 9 ° e p~ . 

deal with a well localized atom we can assume 9 = tt in 
Eq.([9]), which is realized by choosing a proper value of 
e corresponding to the position of atom in the x direc- 
tion. We pointed out that for the experimentally feasible 
parameters of the system under consideration [5TJ] . i.e., 
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fc ~ 10 6 m _1 , m = lOpg, w m /27r = 10MHz, L = 1/im where A 
and for the smallest achievable value of the cavity beam 
waist Wo > 10~ 9 m the LDP is always less than one, i.e., 
77 < 1. In this limit we can keep terms up to first or- 
der in the phonon number n b , and safely truncate the 
summation of Eq. ( 14 1 , 



f J= i(n b ) ~ 1 



-n b . 



(18) 



By substituting Eq.(18) into the Hamiltonian (17 1 one 



can write the interaction of Hamiltonian as 

Hi n t = ~ht;o{b + b^a^a + hg^ waa + + a~(vb 
— hTI 9 » b i >n b aa + _|_ (j-a^ribb 



(19) 



Note that the first and second terms of the above Hamil- 
tonian denote the standard tripartite atom-field-mirror 
coupling which recently has been studied in Refs.[T§] and 
[2"T] . The third term denotes an intensity-dependent cou- 
pling among the three subsystems of atom-field-mirror. 
This type of nonlinear coupling is attributed to the spa- 
tial field-mode structure at the position of the atom. 



0/ = uj c — uji is the cavity detuning with respect 
to the laser, G — g M -\/|(c7 z )| and 7 m is the decay rate of 
the vibrational mode of the MR. The motional quantum 
fluctuation b in (t) satisfies the following relations [55] 



(b in (t)bl(t')) = [{n btth ) + l]5(t-t'), 
(blm n (t')) = {n b , th )S(t-t'), 
(b in (t)b in (t')) = {b\ n (t)b\ n (t')) = 0, 



(21) 



where (n bt th) is the mean number of phonons in the 
absence of optomechanical coupling, determined by the 
temperature of the mechanical bath T, 



{n b fh) 



(22) 



e "b 1 



The only nonvanishing correlation function of the noises 
affecting the atom and the cavity field is (a in (t)a\ n (t')) — 
(F a (t)FKt')) = S(t-t'm- 



III. DYNAMICS OF THE SYSTEM 



To describe the dynamical behavior of the system 
under consideration it is necessary to consider the 
fluctuation-dissipation processes affecting the three sub- 
systems. For this purpose, we first assume the exci- 
tation probability of the single atom to be small. In 
this limit, the dynamics of the atomic polarization can 
be described in terms of the bosonic operators c and 
ct([c,ct] = 1) [24, 52], where the atomic annihilation op- 
erator is defined as c = a~ / yf\{o z )\. This is valid in the 
low atomic excitation limit, i.e., when the atom is ini- 
tially prepared in its ground state [24]. This means that 

the single-atom excitation probability should be much 

1 1 2 

smaller than one. i.e., ^j"^ 2 -C 1, where A a = w a — cui 
is the atomic detuning with respect to the laser and 7 a 
is the decay rate of the excited atomic level. Therefore, 
the bosonization of the atomic operators is valid only if 
go <C A 2 + -fl that is the atom is weakly coupled to the 
cavity. 

The dynamics of the system is fully characterized by 
the following set of nonlinear quantum Langevin equa- 
tions, written in the frame rotating at the input laser 
frequency, 

2 

iA„]c-iG(l- \n b )atf 



A. Linearization of QLEs 



Our aim is to study the conditions under which one can 
efficiently correlate and entangle the atom and the me- 
chanical resonator by means of the common interaction 
with the intracavity optical mode. As shown in Refs. [45] 
and [54] . a straightforward way for achieving stationary 
and robust entanglement in continuous variable optome- 
chanical systems, is to choose an operating point where 
the cavity is intensely driven so that the intracavity field 
is strong, which is realized for high-finesse cavities and 
enough driving power. Therefore, we focus onto the dy- 
namics of the fluctuations around the classical steady 



'[7a 



iA, 



o/J 



i£ a(6 + 6 f ) - iG(l - yn fc )6c + 

o 

- i£o<3^a — iG[(l — rfn b )ac' - 



-[y m + ica, 



state by decomposing each operator in Eqs. (20 1 as the 
sum of its steady-state value and a small fluctuation, e.g., 
a = a s + da, b = f3 s + 5b, and c = c s + 5c. The steady 
state terms of these operators are given by 



b s = 



(uj m - i7 m ) 
G 2 a s 



(20a) 

(20b) 
a f c6 2 ] + 
(20c) 



(i 7a -A a ) 



£ = a. 



iA/ + k 



|G 2 I 2 



( 7a + iA a ) 



(23a) 
(23b) 
(23c) 



where A/ = Ao/ — 2£o-Re(6 s ) denotes the effective op- 
tomechanical detuning and £ = 2^oa s . The other param- 
eters are defined in the appendix. In the linearization 
manner, we also obtain the following linear QLEs for the 



G 



quantum fluctuations of the triple system 
Sc=-[ 7a + iA a ]5c-iG[(l- r ^\b 2 s \)(a s 6tf + b* s Sa) 



^a&ibjtf + b*Sb) 



V2jaF a , (24a) 
da = -[« + iA f]Sa + i£o Sa(b s + b* s ) + a s {Sb + Stf) 

2 2 

- iG[(l - |-|&| 2 )(c 5 .5& + Mc) - yft.C.ffi.^ + K5b) 
2^a in , (24b) 
[7 m + iu> m ]8b + i^ a s (^a t + Sa) 
{l-rj 2 \b s \ 2 )(a s 5J + c*6a)-ri 2 a 3 c* s (b s 6tf +b*5b) 



IG 

,2 



- y{&s(c^a f + a.Se) + 2a s b s c s 6b}\ + y^A^ (24c) 
in terms of the fluctuations of the quadrature operators, 

SX a = ^-={5a + 8a r ),8Y a = -^{5a-Sa^), (25) 
v 2 v 2i 

5X C = -=(5c+6ci),8Y c = -=;(8c-6ci), (26) 
V 2 v 2i 

*g=4=(<y6 + <5& t ),<yp=-7=r(^-<y6 t )- (27) 
v2 v 2i 

The resulting linearized QLEs can be written in the fol- 
lowing compact matrix form 



ii(t) = Au(t) + n(t), 



(28) 



where u(t) = [5q(t),5p(t),5X a (t),5Y a (t),6X c (t),5Y c {t)] T 
is the vector of CV fluctuation operators and 
n(f) = [y^q in (t),y/2^ P in (t),V^X^(t),V2^Y^(t) : 
^/2j a X , c n (t), v /2j a Y t i n (t)] T is the corresponding vector of 
noises. Moreover, the drift matrix A is a 6 x 6 matrix 



The steady state is a zero-mean Gaussian state due to 
the fact that the dynamics of the fluctuations is linearized 
and all noises are Gaussian. As a consequence, it is fully 
characterized by the 6x6 stationary correlation matrix 
(CM) V, with matrix elements 



v _ = (^(00)^(00) + Mj(00K(00)) ^ 



The formal solution of Eq. ( 28 1 yields 



Vij = 



ds ds'M lk {s)M 3 i{s')D kl { s - s'), (31) 
Jo 



where M(t) = exp(At) and D(s - s') is the diffu- 
sion matrix, the matrix of noise correlations, defined as 
D k i(s-s') = (n k (s)ni(s r ) + ni(s')n k (s)) /2. For the noise 
diffusion matrix we have D(s — s') = D6(s — s'), where 
D = diag[7 m (2n b + l),-f m (2n b + 1), k, k, J a , J a ). There- 
fore, Eq.(31) is simplified to 



where 



V= I dojV(uj), 
la 



V(uj) = M{oj)DM{lu) t . 



(32) 



(33) 



When the stability conditions are satisfied (M(oo) = 0), 
one can obtain the following Lyapunov equation 



AV + VA T = -D. 



(34) 



Equation (34) is a linear equation for V and can be 
straightforwardly solved. However, the explicit form of 
V is complicate and will not be reported here. 
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(29) 

where Of and Of denote the real and imaginary parts of 
parameter Oi, respectively. The other matrix elements 
are defined in the appendix. 



B. Stationary quantum fluctuations 

Here, we focus our attention on the stationary proper- 
ties of the system. For this purpose we should consider 
the steady state condition governed by Eq. (28 1. The 



steady state is reached when the system is stable, which 
occurs if and only if all the eigenvalues of the matrix A 
have a negative real part. These stability conditions can 
be obtained, for example, by using the Routh-Hurwitz 
criterion 1 55 . 



IV. ENTANGLEMENT PROPERTIES OF THE 
STEADY-STATE OF THE TRIPARTITE SYSTEM 

In this section we examine the entanglement properties 
of the steady state of the tripartite system under consid- 
eration. For this purpose, we consider the entanglement 
of the three possible bipartite subsystems that can be ob- 
tained by tracing over the remaining degrees of freedom. 
Such bipartite entanglement will be quantified by using 
the logarithmic negativity |56j . 



Epf — max[0, — ln27y ], 



(35) 

r i V2 

where rf = 2" 1/2 £(H P ) - \/^{V bp ) 2 - 4detV &p is 

the lowest symplectic eigenvalue of the partial transpose 
of the 4x4 CM, Vb p , associated with the selected bipar- 
tition, obtained by neglecting the rows and columns of 
the uninteresting mode, 



Vbp - [ c t B i 



(36) 
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FIG. 3. (Color online) Plot of En of the three bipartite sub- 



(atom-mirror), Ej^ (field-atom 



t— rTTT f { 

E N (mirror- 

field)]versus the normalized atomic detuning A a /u m at fixed 
temperature T — 0.4 K, and for two different values of the 
LDP: (a) r\ = 0.04 and (b)?, = 0.08 . The optical cavity 
detuning has been fixed at A; = — uj m , while the other pa- 
rameters are aj m /27r = 10 MHz, Q = 11 X 10 5 , m = 10 pg, 
k = 0.07u m , k = 10 6 m _1 , P c = 800fiW, j a /2TT = 0.04w m , 
and go/2% = 10 3 Hz. 



and S(Vbp) = dctB + detB' - 2detC. 

In Fig. 3 we have plotted the three bipartite logarith- 
mic negativities, EfJ 71 (atom-mirror) , (field-atom), 
and Ex (mirror-field) versus the normalized atomic de- 
tuning A a /oj m at fixed temperature T — 0.4 K, for 
two different values of the LDP [ r; = 0.04, Fig. 3a 
and rj = 0.08, Fig. 3b] and for the experimentally fea- 
sible parameters|50). i.e., a mechanical resonator with 
oscillation frequency uj m /2ir = 10 MHz, quality factor 
Q = 11 x 10 5 , m = 10 pg and an optical cavity with 
length L = 1 /im and damping rate k — 0.07w m driven 
by a laser with kg ~ 10 6 m _1 and power P c — 800 fiW. The 
atom damping constant has been taken 7 a /27r = 0.04w m 
with coupling constant go/2ir = 10 3 Hz. 

The optical cavity detuning has been fixed at Af = 
—uj m which turns out be the most convenient choice. 
As seen, by increasing LDP, the two bipartite entangle- 
ment of Efp and E f N a increase and the bipartite entan- 
glement of E^ decreases overall. The reason is that, 
by increasing the LPD the tripartite atom-field-mirror 
coupling rate increases compared to the coupling rate of 



the bipartite field-mirror subsystem, or equivalently the 
parameter G/£ is increased. This result reveals that by 
changing the LDP one can control the tripartite coupling 
amplitude or even goes through the regime in which the 
tripartite system reduces to an effective bipartite subsys- 
tem. However, the three logarithmic negativities do not 
behave in the same way and the entanglement sharing is 
evident. In particular, the entanglement of interest, i.e., 
increases at the expense of the mirror-field entan- 



77 am 



glement, while E^ 1 remains always non- negligible. 




-1.0 -0.5 0.0 
A a /(O m 



1.0 



0.8 
0.6 
feJ 0.4 
0.2 
0.0 





, , , 1 , , , 


11 


11 


(b) : 




f 1 


nam 

-N 









































-2.0 -1.5 -1.0 -0.5 0.0 
AJ<o m 



0.5 



1.0 



FIG. 4. (Color online) Plot of En of the three bipar- 
tite subsystems [Eff 1 (atom-mirror), E^ (field-atom), E™ 
(mirror-field)]versus the normalized atomic detuning A a /tj m 
for a fixed value of the LDP, rj — 0.04, and for two different 
temperatures: (a)T = 1.2K and (b) T = 3K. The optical 
cavity detuning has been fixed at Af = — uj m and the other 
parameters are as in Fig. 3. 

Fig. 4 shows the logarithmic-negativity of the three bi- 
partite subsystems versus the normalized atomic detun- 
ing A a /u) m for a fixed value of the LDP, 77 = 0.04, and 
for two different temperatures: T = 1.2 K (Fig. 4a) and 
T = 3 K (Fig. 4b). The optical cavity detuning has been 
again fixed at Af — —uj m . As expected, the three kinds 
of bipartite entanglement decrease by increasing temper- 
ature, but the atom-mirror and field-atom entanglement 
show highly temperature robustness. However, the field- 
mirror entanglement shows extremely fragile entangle- 
ment robustness versus temperature and its logarithmic- 
negativity falls down to zero at T = 3K. 

Generally, the scheme is able to generate apprecia- 
ble entanglement between the atom and the MR, es- 



pecially by sharing from the mirror-field entanglement. 
Similar bipartite entanglement behavior can be observed 
in other similar tripartite systems, such as the atom- 
field-mirror scheme proposed in Ref. |24j , the microwave- 
optical-mirror system of Ref. |59| and the two cavity op- 
tomechanical setup of Ref. [60] . 
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FIG. 5. (Color online) Plot of En of the three bipartite sub- 
systems [Efj 1 (atom- mirror), E^ (field-atom), E^ (mirror- 
field)] versus the LDP, r), for fixed temperature T = 0.4K. 
The optical cavity and atomic detuning are Af = —co m and 
A a = uj m , respectively. The other parameters are as in Fig. 
3. 

The effect of the LDP is also illustrated in Fig. 5, where 
we have plotted the logarithmic-negativity as a function 
of 1] at fixed optical cavity detuning Af = —u> m and at 
atomic detuning A a — uj m . ft is clear that by increasing 
the LDP, the field-atom entanglement increases when the 
entanglement between the intracavity mode and the MR 
is drastically suppressed. We also observe that the atom- 
mirror logarithmic-negativity EfJ™ plateaus in a case as 
r\ increases. 

? fa and Eff 1 



Fig.6 shows E J N 



versus A a /uj m and ^ a /^n 



fa, 
N 



and E% rn 



are max- 



for r] — 0.04. As is clear, the E' 
imized around sideband A a ~ uj m . By increasing the 
atomic spontaneous emission rate j a as we expected, 
both logarithmic-negativities decrease drastically. 

The entanglement features of the tripartite system at 
the steady state can be observed by experimentally mea- 
suring the corresponding CM. This can be done by com- 
bining existing experimental techniques. By homodyn- 
ing the cavity output one can measure the cavity field 
quadratures. Ref. [3] has proposed a scheme to measure 
mechanical position and momentum of the MR, in which 
by adjusting the detuning and bandwidth of an addi- 
tional adjacent cavity, both position and momentum of 
the mirror can be measured by homodyning the output 
of the second cavity. Moreover, by adopting the same 
scheme of Ref. [61] . the atomic polarization quadratures 
X a and Yq can be measured, i.e., by making a Stokes pa- 
rameter measurement of a laser beam, shined transversal 
to the cavity and to the cell and off-resonantly tuned 
to another atomic transition. Very recently, Ref. [52] has 
demonstrated the proof of principle of the use of a Bosc- 
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FIG. 6. (Color online) Density plot of En of the bipartite 
subsystems: (&)E^ and (b) _E^, m versus A a /oj m and 7 a /oj m 
for r\ = 0.04 and for T — 0.4K. The optical cavity detuning 
has been fixed at A/ = — u) m ■ The other parameters are as in 
Fig. 3. 



Einstein condensate(BEC) as a diagnostic tool to deter- 
mine the elusive mirror-light entanglement in a hybrid 
optomechanical device. In such a case, one dose find a 
working point such that the mirror-light entanglement is 
reproduced by the BEC- light quantum correlations. 



V. NORMAL-MODE SPLITTING IN THE 
DISPLACEMENT SPECTRUM OF THE MR 

In this section, we show that the atom-field-mirror cou- 
pling leads to the splitting of the normal mode into three 



9 




modes [Normal Mode Splitting(NMS)]. The optomechan- 
ical NMS however involves driving four parametrically 
coupled nondegenerate modes out of equilibrium. The 
NMS does not appear in the steady state spectra but 
rather manifests itself in the fluctuation spectra of the 
mirror displacement. To study the NMS in our system 
we need to find out the displacement spectrum of mirror 
as: 



Sg(u) 



2tt 



dQe-^+^iSq^SqiQ) + Sq((l)Sq(u)) 



(37) 



where Vu(uj,fl) = l/2(Sq(uj)5q(il) + 5q(Q)5q(ui)) is an 
element of CM which is given by Eq.(33). Unfortunately, 



the analytical form of the displacement spectrum of the 
mirror is too complicate to put a clear physical interpre- 
tation on it. Thus, in the following, we give and analyze 
the results obtained by numerical calculations. 




FIG. 7. (Color online) Normalized plot of the displacement 
spectrum S q (cj) versus ui/ui m at fixed temperatures T = 0.4K 
and for two different values of the LDP: r\ = 0.016(red-line), 
77 = 0.04(blue-line). The optical cavity and atomic detuning 
have been fixed at A/ = ui m and A a — uj m , respectively. The 
other parameters are as in Fig. 3. 

Fig. 7 shows the displacement spectrum of the MR as 
a function of the normalized frequency Lo/uj m at Af = 
u) m , A a = u> m and for two different values of the LDP: 
r\ = 0.016, r\ = 0.04. For the small values of the LDP, we 
observe the usual normal-mode splitting into two modes 
with central peaks at the sidebands ui = ±cJ m . This 
figure shows a highly symmetric structure with respect 
to ui = 0. As is seen, by increasing the LDP the normal 
mode splits up into three modes. 

A more clear illustration of the three-mode splitting 
is shown in Fig. 8. This figure shows the displacement 
spectrum of the MR versus the normalized frequency 
u;/u; m and atomic detuning A a /uj m at Af = u) m . The 
three-mode splitting manifests itself mainly at A a ~ uj. 
By going through the region far from A a ~ u), three- 
mode splitting merges into two-mode splitting around 
A a ~ 0.75a; m and A a ~ 1.35w m . The NMS is associated 
with the mixing among the vibrational mode of the MR, 
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FIG. 8. (Color online) Density plot of the displacement spec- 
trum S q (uj) versus uu/uim and A a /uj m for T = 0.4K and 
77 = 0.04. The optical cavity detuning has been fixed at 
Af = uj m . The other parameters are as in Fig. 3. 



the fluctuations of the cavity field around the steady state 
and the fluctuations of the atomic mode. The origin of 
the fluctuations of the cavity field is the beat of the pump 
photons with the photons scattered from the atom. For 
not so large values of the LDP (small nonlinearity) the 
field-atom coupling is much smaller than the field-mirror 
coupling. Therefore, the system simply reduces to the 
case of two mode coupling, i.e., coupling between the 
mechanical mode and the photon fluctuations [34] . When 
the LDP is large enough the mechanical mode, the pho- 
ton mode, and the atomic mode forms a system of three 
coupled oscillators. The occurrence of splitting of the 
normal mode into three modes has been analyzed re- 
cently in another tripartite system, i.e., a cavity quantum 
optomechanical system of ultracold atoms in an optical 
lattice 34 . Furthermore, similar three coupled oscilla- 
tor experimental results where two coupled cavities, each 
containing three identical quantum wells [63] and one mi- 
crocavity containing two quantum wells [64] have been re- 
ported. 

It should be pointed out that to observe the NMS, the 
energy exchange between the three modes should take 
place on a time scale faster than the decoherence of each 
mode. The normal mode splitting into three modes due 
to local increasing of the LDP has also been reported in 
Ref.[BS], where the authors have shown that the NMS 
can be observed only if the coupling between the atoms 
and the cavity is strong enough. This strong coupling 
can be achieved by increasing the atom numbers. One 
experimental limitation could be spontaneous emission 
which leads to momentum diffusion and hence heating of 
the atomic sample [66]. In our model, we don't encounter 
such a limitation and three-mode splitting is approached 
by proper choosing of the LDP. 
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VI. CONCLUSION 



the implementation of quantum teleportation and/or the 
photon blockade process to prevent two or more photons 
from entering the cavity at the same time. 



In this work, we have proposed a theoretical scheme 
for the realization of tripartite intensity-dependent cou- 
pling among a single mode of a Fabry-Perot cavity with 
an oscillating mirror, a single two-level atom inside it, 
and a vibrational mode of the oscillating mirror. We 
have shown that in the presence of Gaussian standing- 
wave of the optical cavity mode, a type of tripartite be- 
tween atom-mirror-field coupling can be manifested. To 
describe such interaction we then have found the general 
form of the corresponding nonlinear Hamiltonian. We 
have restricted our investigation to first vibrational side- 
band j = 1 and have studied its dynamics by adopting 
a QLE treatment. We have focused our attention on 
the steady state of the system and in particular, on the 
stationary quantum fluctuations of the system by solv- 
ing the linearized dynamics around the classical steady 
state. We have seen that, in an experimentally acces- 
sible parameter regime, the steady state of the system 
shows both the tripartite and the bipartite CV entan- 
glement. We have shown that the LDP(as a measure of 
the strength of nonlinearity)can extremely modifies both 
the tripartite and the bipartite CV entanglement in the 
system. In particular, by increasing the LDP, one can 
see that the field-atom and atom-mirror entanglement 
increase at the expense of optical-mechanical entangle- 
ment. The intracavity mode is able to mediate for the 
realization of a robust stationary (i.e., persistent) entan- 
glement between the MR mode and the single two-level 
atom, which could be extremely useful in quantum in- 
formation/quantum computer networks in which the MR 
modes are used for quantum communications [67, 68 , and 
the atom is used as a qubit(e.g., solid-state qubits). Fur- 
thermore, we have analyzed the occurrence of the NMS in 
the displacement spectrum of the oscillating mirror. As 
we have shown, for a small value of the LDP, the usual 
normal-mode splitting into two modes with central peaks 
at the sidebands u> — ±oj m is observed and by increasing 
the LDP the normal mode splits up into three modes. We 
have shown that, when the LDP is large enough the me- 
chanical mode, the photon mode, and the atomic mode 
forms a system of three coupled oscillators. The realiza- 
tion of such a scheme will also open new opportunities for 
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Appendix: Definition of the elements of the drift 



matrix of Eq.(29) 



In the drift matrix of Eq. ( 29 ) we have defined 



Tim = 7m + M3 , r 2m = 7 m + M\ , 
film = w m + M$ , Q,2rn = &m + M± i 



Gi=G 



c ;(i~ v i\ bs \ 2 ) + r Lc s b 2 s 



G 2 = G&:(1-||6 S | 2 ), 



G 3 = G 
M x = -Ga 
M 2 = Ga, 
M 3 = Gc s 



c:{i-v 2 \bs\ 2 )-\c s bi 



(l-V 2 \bs\ 2 )+ V Y bf 

(l-r? 2 |6 s | 2 )-^r 



(l-V 2 \b s \ 2 ) + \bl 



M 3 = Gc s 
Mi = -Grfa. 
M 5 = -Gifa. 



6*c* + c s b s - c* s b s 



b* s c* + c s b s + c*b s 



N 1 = b s (i- r -L\b s \ 2 ), 



N 2 = - 
N 3 =a s 



(l-V 2 \b s \ 2 ) + \bi 

(i-v 2 \b s \ 2 )-^bi 
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